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SUMMARY 



Results are presented of drag tests of six bodies of 
revolution with systematically varying shapes and with a 
fineness ratio of 5. The forms were derived from source- 
sink distributions, and formulas are presented for the 
calculation of the pressure distributions of the forms. 
The tests were made" in the N.A.C.A. variable-density tun- 
nel over a range of values of Reynolds Number from about 
1,500,000 to 25,000,000. The results show that the bodies 
with the sharper noses and tails have the lowest drag co- 
efficients, even when the drag coefficients are based on 
the two-thirds power of the volume. The data show the 
most important single characteristic of the body form to 
bo the tail angle, which must be fine to obtain low drag. 



INTRODUCTION 



The National Advisory Committee for Aeronautics is 
conducting an investigation of fuselage drag at large 
values of the Reynolds Number. The first phase of the in- 
vestigation is the study of streamline bodies of revolu- 
tion to obtain data on basic forms that can be modified to 
obtain fuselages of practical shape for particular types 
of airplanes. 

Tests were made of six streamline bodies of revolu- 
tion with systematically varying shapes and a fineness 
ratio of 5. The tests were made in the fall of 1935 in 
the N.A.C.A. variable-density wind tunnel. The results 
cover a range of Reynolds Numbers from about 1,500,000 to 
25 ,000 ,000. 
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DERIVATION OF STREAMLINE FORMS 



The streamline forms were derived from systematically 
varied source-sink distributions (references 1 and 2) cov- 
ering a range of shapes in which it was thought that the 
lowest drag forms for the selected fineness ratio of 5 
would "be found. The type of source-sink distribution used 
is shown in figure 1. The abscissas of the sources are 
denoted by £ , the intensity of the sources per unit 
length by f(£) in which positive values of f(|) denote 
sources and negative values, sinks. The condition that 
makes the stream from the sources self-contained is ex- 
pressed by the equation 



1 f(t> <U = o (D 



The stream function is taken as the quantity of fluid 
flowing inside a circle perpendicular to the axis and hav- 
ing its center on the axis (fig* 1). In the computation 
of the stream function \j/ due to the sources, the amount 

of fluid delivered by the sources upstream from the circle 
is deducted from the quantity passing through the circle. 
The stream function of the sources is then given by the 
formula 

\ - - § {"fit) (l + i-=J.) dt (2) 

0 

where r = J (x - if + y 2 (reference 2). 

To this stream function must be added that due to the 
superimposed parallel flow 

\J/ S = V tr y a (3) 

The condition defining th-c surface of the body about which 
the flow takes place is 



\|^ + v[/ a = * = 0 (4) 
The source distributions used are shown in figure 
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2(a,"b) and the defining constants are given in table I. 
The integration of equation (2) for the type of source 
di stri tmtion used and the condition 



,f 



f(| ) a| = l J fit ) d£ = - i 



0 g 

results in the following equation: 

1 r 
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The formulas for the speeds of flow due to the source 
distribution normal to the axis? (v) and ijarallel to the 



axis (u) are as follows: 



1 j' 2 ^xe - x 3 - y 2 + 2 + y 2 ^ + 2 /f - x e - x\ 
4-rry I of V r Q i' / f - e \ r T » r„ / 



(f - e) 



,^xf_ -_x a __-_y 2 xe_-_x^_-_y 2 \ + 2g 

^ i> '* r / (h - g)(l - 



h - x g - x\ 



^ r h 



,2 2 3 2 

xh - x - y xg - x - - 



r„ ; (h - g)(l - e) v 



e) 



g 



/"2x - x 2 - y 3 - 1 xh - x 2 - y 2 - h + x N ^ 



(1 - g)(l - h) 



)} 



(6) 
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It is evident that a group of "body surfaces of differ- 
ent fineness ratios can "be obtained from one source dis- 
tribution "by choosing different values for the superimposed 
parallel flow. The maximum diameter of the bodies to be 
■derived was arbitrarily chosen as one-fifth of the length 
of the source distribution. The speed 7 of the superim- 
posed parallel flow to satisfy this condition was deter- 
mined as follows: 

The station x of the maximum diameter of the body 
was found by determining from equation (6) the value of x 
at which the velocity v is zero at a value of y corre- 
sponding to the selected maximum radius of the body. The 
stream function was then computed for these coordi- 

nates by means of equation (5). The required value of V 
was then readily determined from eqtiations (3) and (4). 

The coordinates of the surfaces of the bodies given 
in table II were obtained by computing a number of values 
of and \|/ a corresponding to selected values of x 

and y, and graphically solving equation (4). 

This graphical solution (reference 2) consists of 
plotting curves of \J/ X for the condition x = constant 
against y« The intersections of these curves with that 
of equation (3) determine the coordinates x and y of 
points on the surface of the body. 

In general, the nose and tail points of the body will 
not fall at x = 0 and x = 1. For the range of source 
distributions used, however, the tail points were not ap- 
preciably displaced from the point x = 1 and are consid- 
ered to occur at that point. The nose points, however, 
are appreciably upstream from the point x = 0 and are 
not readily determined with sufficient accuracy by the 
foregoing methods. These points were determined by finding 
the point on the x axis where the velocity from the 
source distribution is equal and opposite to that of the 
superimposed flow. This point was determined from a sim- 
plified form of equation (7) derived to hold on the x 
axis ii'i the neighborhood of x = 0. The fineness ratios 
of the forms obtained are slightly greater than 5 because 
of this extension of the nose ahead cf the point x = 0# 

The ordinates of the ten forms computed are given in 
table II and the outlines of the forms are shown in figure 
3. Each form is designated by a number of three digits 
which indicate in a general way the nose shape, the tail 
fullness, and the tail angle. 
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The theoretical pressure di st r i "but i on over the surface 
of the "body or in the field about the body can "be obtained 
from Bernoulli ! s equation and the velocities v and u 
computed from equations (6) and (7). The total velocity 
at any point is found "by adding vectorially these veloci- 
ties to the velocity V of the superimposed parallel flow. 
The theoretical pressure distribution about form 111 has 
been computed and is shown in figure 4. It can be seen 
that the points of, discontinuity of the source-distribution 
curve have no marked direct effect on the fairness of the 
pres sur e-di stributicn curve- 



MODELS AND TESTS 



Aluminum-alloy models were made of forms 111, 221, 
222, and 332. The models were 8 inches in maximum diame- 
ter and were carefully finished to prevent drag increases 
due to surface roughness (references 3 and 4). Each model 
was built in two sections, the division being made at the 
maximum diameter to allow the various nose and tail por- 
tions to be combined to form models of different shape. 
In this way models approximating the shapes of any of the 
remaining six computed forms could be obtained if tests of 
these forms were found to be desirable as a result of the 
tests of the original four models. Tests were actually 
obtained in this way of forms 121 and 211. The other four 
forms were not tested because the results indicated that 
they would have less favorable drag characteristics than 
the best tested forms* 

The models were tested at zero yaw and pitch in the 
N.A.C.A. variable-density wind tunnel (reference 5 ) at 
six values of the Reynolds Number, based on the model 
length, from .about 1,500,000 to 25,000,000. The test 
methods and corrections applied to the results are de- 
scribed in references 3 and 4. The precision of the tests 
was as described in reference 4 except that the uncertainty 
of the balance calibration mentioned therein had been elim- 
inated by the use of an improved drag balance. 

RESULTS A"TD DISCUSSION 

Drag coefficients for the six models tested based 
both on the cross-sectional area C^ and the two-thirds 

^A 
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power of the volume C-n are shown plotted against Reynold 
Number in figure 5. The "bodies with the sharper noses and 
tails have the lowest drag coefficients even when the co- 
efficients are "based on the two-thirds power of the vol- 
ume. The data show the most important single characteris- 
tic of the "body form to he the tail angle, which must he 
fine to obtain low drag. 

Values of the drag coefficient corrected to an effec- 
tive value of the Reynolds Number of 66,000,000 from the 
test value of the Reynolds Number of 25,000,000 are tabu- 
lated in table III. The correction (reference 6) allows 
for the decrease in skin-friction drag at the effective 
value of the Reynolds Number below that for the test value 
of the Reynolds Number. The correction is made by multi- 
plying- the tost values of the drag coefficients by a fac- 
tor, which is taken as 0,875 for a test value of the 
Reynolds Number of 25,000,000. The corrected values of 
the drag coefficient are believed to be more nearly appli- 
cable to flight at the effective Reynolds Number than the 
test values at the test Reynolds Number. 

The fuselage size for some airplanes, particularly 
small ones, is largely determined by the cro ss- sec tional 
area required by the selected seating arrangement. In 
such cases a fuselage shape such as 111 may save an appre- 
ciable part of the fuselage drag as compared with a more 
conventional shape such as 222, In the case of other air- 
planes, such as large transports, the longitudinal distri- 
bution of fuselage volume becomes of considerable impor- 
tance, and the fuselage size may be determined by consid- 
erations other than that of maximum cro s s- s ect i onal area. 
Under these circumstances a careful analysis is required 
to select the proper form to give lowest drag. Such an 
analysis is hampered by the lack of knowledge of several 
factors including the variation of drag with fineness rati 
a study of which constitues the next part of the N.A.C.A. 
investigation of fuselage drag. 



Langlcy Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , August 5, 1937. 
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TABLE I 

CONST A.: I DEFINING THE SOURCE-SINK DISTRIBUTIONS USED 



Form 


c 


f 


r, 


h 


111 


0.050000 


0. 333333 


0.400000 


0.700000 


- o *: 
J.. -J. 


.05 0000 


. 3 33333 


.5 25 000 


.775000 


122 


.05 0000 


. 333 333 


.450 000 


.925000 


211 


.012500 


. 25 0000 


.400000 


.700000 


221 


.012500 


. 25 0000 


.525000 


. 775000 


222 


.012500 


. 250000 


.450000 


.925000 


232 


.012500 


. 25 0000 


.550000 


.950000 


321 


,006250 


. 175 000 


.525000 


.775000 


322 


.006250 


.175000 


.450000 


.9 25000 


332 


.006250 


.175000 


.550000 


.950000 
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TABLE III 
DRAG- COEFFICIENTS AT AN EFFECTIVE 
REYNOLDS NUMBER OF 66,000,000 
(Reynolds Nunher is "cased on fuselage length) 



Form 




\ 


111 


0.0179 


0.0401 


121 


.0176 


.0405 


211 


.0177 


.0415 


2 21 


.0178 


.04 37 


22?. 


.0186 


.0472 


332 


.0193 


.0508 
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Figure 5.- Drag ooaffioienta of K.A.C.A. fuselage foma. 



